function y_new = rk4(y0, dt, param)
    % 使用四阶Runge-Kutta方法进行单步更新
    % model: 模型函数句柄，定义了系统状态变化
    % y0: 初始状态
    % dt: 步长

    % 计算四阶RK法的中间结果
    k1 = resistive(y0, param);                         % 第一阶段（k1）
    k2 = resistive(y0 + 0.5*k1, param);                % 第二阶段（k2）
    k3 = resistive(y0 + 0.5*k2, param);                % 第三阶段（k3）
    k4 = resistive(y0 + k3, param);                    % 第四阶段（k4）

    % 更新状态
    y_new = y0 + dt/6*(k1 + 2*k2 + 2*k3 + k4)';       % 使用RK4公式更新状态
end
